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ABSTRACT 

Context. The baryonic Tully-Fisher relation (BTF) is a fundamental relation between baryonic mass and maximum rotation velocity. 
It can be used to estimate distances, as well as to constrain the properties of dark matter and its relation with the visible matter. 
Aims. In this paper, we explore if extremely low-mass dwarf galaxies follow the same BTF relation as high-mass galaxies. We quantify 
the scatter in the BTF relation and use this to constrain the allowed elongations of dark matter halo potentials. 

Methods. We obtained H i synthesis data of 1 1 dwarf galaxies and derive several independent estimates for the maximum rotation 
velocity. 

Results. Constructing a BTF relation using data from the literature for the high-mass end, and galaxies with detected rotation from 
our sample for the low-mass end results in a BTF with a scatter of 0.33 mag. 

Conclusions. This scatter constrains the ellipticities of the potentials in the plane of the disks of the galaxies to an upper limit of 
0-0.06, indicating that dwarf galaxies are at most only mildly tri-axial. Our results indicate that the BTF relation is a fundamental 
relation which all rotationally dominated galaxies seem to follow. 

Key words, dark matter - galaxies: kinematics and dynamics - galaxies: dwarf - galaxies: fundamental parameters 



1. Introduction 

■ The Tully-Fisher (TF) relation dTullv & FisheJll977l) . the rela- 
tion between the luminosity and rotation velocity of a galaxy, 
has been extensively use d to estimate extra galactic distances 
(e.gJPierce & Tullvjl98&IWiliick et al.ll997t lSakai et al.l2000t 
iTully & Pierce 2000; Sp ringob et alJl2007l) . In the usual "classi- 
cal" interpretation, luminosity is a proxy for the stellar mass, 
which in turn depends on the total (visible and dark) mass 
and through it on the rotation velocity. The slope and zero 
point of this classical TF relation do not depend on the cen- 
tral surface brightness of galaxies (IZwaan et al i ri995h . though 
for v ery low mass dwarf galaxies, the slope tend s to steepen 

■ (e.g., iMatthews etail Il998t iMcGaugh etafl |2000|) . Low-mass 
dwarf galaxies are apparently underluminous given their rota- 
tion velocity and therefore fall below the TF relation as de- 
fined by the high mass galaxies. A single linear relation can 
be restored if one replaces the luminosity (or stellar mass) with 
the baryonic disk mass, thu s including the gas mass as well 
dMcGaugh et al.ll999ll2000h . This relation is called the baryonic 
Tully-Fisher (BTF) relation an d has been studied by many au- 
thors i n the last few years (e ,g., Bell & de Jondl200ltlVer heiien 
200 1 [ iGurovich et al.ll2064jMcGaughll2004 120051; beha et al.1 
20061; iDe Riicke et alj 120071; iNoordermeer & Verheijenl 12007; 
Star k et alJl2009h . 

The existence of a (baryonic) Tully-Fisher relation places 
severe constraints on galaxy formation and evolution the o- 
ries (cf.lEisenstein & Loeb1ll996tlMcGaugh & de Bloklll998allbl ; 
iMo et alJl998tlSteinmetz & Navarro! 1999l;lBlanton et al.l2 008). 
For example, iFranx & de Zeeuwl (1 19921) note that the observed 
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scatter in the TF relation places upper limits on the elongation of 
dark matter halos. They find that that the ellipticity of the poten- 
tial in the plane of the disk is most li kely between and ~ 0.06 . 
This is in good agreement with what Trach ternach et al.1 {2008) 
found observationally for a sa mple of 1 8 dwarf and spiral galax- 
ies fr om the THINGS survey dWalter et alj|2008t Ide Blok et af] 
2008), but in disagreement with results from CDM simulations 
which predict more elong ated potentials (e.g jFrenk et al.l Fl988; 
lHava shi & Navarro 200 6J). The sm aller scatter in the observed 
BTF relation (e.g., [McGaugh 2005), means that it can be used to 
put similar constrai nts on lower ma ss galaxies as well. This was 
recently shown bv lStark et al.l d2009l) who studied the BTF for a 
large sample of dwarf galaxies with resolved H i rotation curves, 
and found a linear BTF relation with small scatter. 

Th e low mass end o f the BTF re lation has also been stud- 
ied bv lGeha et a l. (2006) and lKovad d2007l) . They measure line 
widths and after correcting these for broadening due to turbu- 
lent motion of the H i, they find that the extreme dwarf galax- 
ies of their sample follow the same BTF relation as the high 
mass galaxies, albeit with a larger scatter. This increas ed scat- 
ter is m ost likely a result of their use of the W20 profile. iBroeilsl 
d 1992b and I Verheijenl d 19971) already pointed out that using the 
maximum rotation velocity from a resolved rotation curve sig- 
nificantly decreases the scatter as compared to using line width 
measurements. 

In this paper, we attempt to determine several independent 
estimates for the maximum rotation velocity ( V max ) for a sample 
of extremely low-mass dwarf galaxies. Doing this, we can check 
whether these galaxies follow the same tight correlation between 
baryonic disk mass and rotation velocity as th eir high mass coun - 
terparts. As a reference, we use the sample of McGaugh (2005), 
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for which well-determined estimates for V max exist based on 
analysis of well-resolved rotation curves. For an in-depth de- 
scription of the comparison sample, the reader is referred to 
McGaugh (2005) and refer ences therein. 

The paper is organized as follows: in Section 2 we describe 
the observations and the data reduction, which is followed by a 
description of the different methods in estimating V max in Section 
3. We comment the individual galaxies in Section 4, and present 
further analysis and our results in Section 5. We summarize our 
results and give our conclusions in Section 6. The Appendix 
contains moment maps, position-velocity diagrams and channel 
maps of our sample galaxies. 

2. The data 

Our sample w a s sel ected from the larger sample of 
ISchombert et al.l dl997[) . which is one of the largest sam- 
ples of extreme field dwarf g alaxies for which both lin e width 
measurements and H i masses (Eder & Schombert 2000), as well 
as optical photometry dPildis et al.l 1 19971) exist. The galaxies 
were chosen to be relatively nearby (Vhei < 1400 kms -1 ), have 
suitable optical inclinations for potential derivation of their 
rotation curves (45 < / opt < 75), and to have V- and /-band 
photometry available. 

2.1. Observations 

Observations were carried out at the Westerbork Synthesis Radio 
Telescope (WSRT) in maxi-short configuration, in the 21 -cm 
line of neutral hydrogen. We use both polarizations and sample 
1024 channels with a bandwidth of 10 MHz or 20 MHz (cor- 
responding to velocity resolutions of 2.10 or 4.12 kms -1 ). The 
sample integration time was set to 60 seconds. Further observa- 
tional details are summarized in Table Q] 

2.2. Data reduction 

The calibration and data reduction of the data is performed us- 
ing standard routines in MIRIAE0. The data are calibrated us- 
ing one of the standard primary calibrators used at the WSRT 
(3C48, 3C147, 3C286, J2052+362). The primary calibrators are 
also used for the bandpass and gain corrections. 

The line-free channels (i.e., the channels containing only 
continuum emission) were used to create a continuum image. 
This image was used to self-calibrate the data. After the self- 
calibration and continuum subtractio n, image cube s were created 
using the robust weighting scheme (Briggs 1995}) with a robust 
parameter of zero for all galaxies. These image cubes were then 
cleaned down to a level of ~ lc For a list of the beam sizes 
and noise levels, see Table Q] We additionally created Hanning 
smoothed data cubes, which were used for the creation of the H i 
profiles and the derivation of the velocity widths. 

Using GlPS^fl we created zeroth, first and second moment 
maps of all galaxies. In order to isolate significant signal, we 
smoothed the data cubes to half the original resolution and only 
retained pixels with values > 2.5 times the (smoothed) noise 
value. Spurious pixels were blotted by hand. This smoothed and 
blotted data cube was used as a mask for the original data cube. 
We determined the number of channels with significant emission 

1 Multichannel Image Reconstruction, Image Analysis and Display 
dSault et al.|[l995l) 

2 GIPSY, the Gr oningen Image Processing SYstem 
dvan der Hulst et al J 19921) 



contributing to each unblanked pixel in the unsmoothed moment 
maps and created a map containing the signal-to-noise (S/N) of 
each pixel in the total intensity (zeroth moment) map. Using this 
map, we determined the average pixel value in the zeroth mo- 
ment map corresponding to a S/N = 3, and clip all moment 
maps using this flux limit. Channel maps and moment maps are 
presented in the Appendix. 

3. Estimating the maximum rotation velocity 

In order to construct a (baryonic) Tully-Fisher relation, one 
needs to estimate the maximum (outer) rotation velocity (V max ) 
of a galaxy. There are several ways in which this can be done, 
which are described below in increasing order of preference. 

3.1. Hi velocity profile 

The simplest way is to use the width W of the global H i veloc- 
ity profile, usually measured at the 20 (50) percent level of the 
maximum intensity, and denoted as W20 (W50). The advantage in 
using the profile width is that it is easy to measure and can be de- 
rived using low resolution data. The drawback is that one cannot 
distinguish between rotation and turbulence. This uncertainty 
matters little for large, fast rotating spiral galaxies. However, 
as one goes to smaller and more slowly rotating galaxies, tur- 
bulent motions will star t to contribute significantly to the total 
width of H 1 profiles (cf. lVerheiienll997l) . In our analysis, we use 
the velocity profiles from the Hanning smoothed data cubes and 
correct them for instrumental broadening and turbulent motion. 
The corrections applied are addressed and discussed more fully 
in Section l5T2l 

3.2. Major axis position-velocity diagram 

An alternative way to derive V max is to make use of the major 
axis position-velocity (pV) diagram. For a galaxy with a flat ro- 
tation curve, Gaussian fits to the outer galaxy radii in the major 
axis /?V-diagram (see, e.g., Fig. lA.ll in the Appendix) yield an 
estimate for the (projected) amplitude of the rotation. The maxi- 
mum rotation velocity is then calculated as half the difference in 
velocity between approaching and receding sides, corrected for 
inclination: 

^receding ^approaching 



The /?V-diagram is potentially able to better constrain the the 
maximum rotat ion curve velocity than W20 and W50 values 
dVerheijenll200ll) . 

3.3. Tilted-ring models 

A third method for the derivation of V max is a tilted-ring model, 
in which the kinematics of a galaxy are described using a set 
of concentric rings. Each of these rings can have its own center 
position (xo,yo), systemic velocity V sys , rotation velocity V rot , 
inclination i, and position angle PA. The pre-requisite for this 
method is a resolved velocity field showing signs of rotation. 

For the galaxies in our sample that meet this criterion, we 
derive rotation curves using the GIPSY task rotcur. Assuming 
that the gas moves on circular orbits, the line-of-sight velocity 
can be expressed as: 

V(x, y) = V sys + V mt (r) sin© cos(0). (2) 



Table 1. Properties of the data. 
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57.7 


12.0 
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1.1 
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2.10 


D640-13 


10-11-2004 


12h 


70.8 


11.5 
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2.10 
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12h 
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0.6 
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" The data for D575-1 are combined from two individual observations 

* The data for D575-2 are combined from three individual observations 
Notes: (1): galaxy identifier; (2): date of observations; (3): length of observations; (4, 5) major and minor axis diameter of the robust weighted 
beam in arcsec; (6): position angle of the beam (in degrees), measured counter-clockwise from the north; (7): noise per channel in mjy beam -1 ; 
(8): pixel size in arcsec; (9): channel width in km s -1 



Here, 6 is the azimuthal distance from the major axis in the plane 
of the galaxy and is related to the position angle PA of the galaxy 
as measured in the plane of the sky by 

-O - x Q ) sin(PA) + (y - y ) cos(PA) 

cos(#) = (3) 

r 

and 

sin(6/) = -(* ~ x ») cos(PA) - (y - y ) sin(PA) 
r cos(;) 

The PA is measured counter-clockwise from the north to the ma- 
jor axis of the receding side of the galaxy. As positions along the 
major axis of a galaxy carry more rotational information than po- 
sitions near the minor axis, we weight the individual data points 
by|cos(0)|. 

The derivation of a rotation curve is generally an iterative 
process involving the consecutive determination of the various 
tilted-ring parameters. Following is a general description of the 
procedure applied. In order to get good initial estimates for i, 
PA, and the center position, we fit isophotes at varying intensity 
levels to the H i total intensity maps, taking care that the results 
are not affected by small-scale structures. As initial estimates for 
V sys and V r ot, we use the central velocities of the W50 profile and 
jWso/ sin(0, respectively. 

In the first tilted-ring fit, we determine the systemic velocity 
by keeping all parameters except V sys and V rot fixed. In a sec- 
ond run, we derive the position of the dynamical center leaving 
by leaving (only) the central position of the rings and V rot un- 
constrained. In a third fit, we determine the position angle by 
leaving only PA, i and V lot as free parameters. As in most cases 
we sample the rotation curves with relatively few (~ 10) tilted- 
rings, we approximate the fitted values either by a constant or 
linearly changing PA. Once the PA is modeled in such a way, 
we make another fit with only i (and V rot ) as free parameters and 
derive the inclination of the galaxy. We do not model any radial 
trends for the inclination. In a last run, we determine the rotation 
curve by keeping all parameters except V rot fixed at their best 
determined values. The maximum velocity of the rotation curve 
derived in this way will be referred to as V v f hereafter. 



4. Comments on individual galaxies 

In this section, we present our results for the individual galaxies. 
Unless mentioned otherwise, the distances given in the following 
sections are calculated by correcting the systemic velocities of 
the galaxi es for Virgocentric infall (following the formalism pre- 
sented in Mo uld et ai]|2000l) . and assuming a Hubble flow with 
a Hubble constant of Ho = 75kms -1 Mpc -1 . Given the differ- 
ent distance outcomes of different flow models, we use half the 
difference between the minimum and maximum flow-corrected 
distances reported in NED as an indication of the uncertainty in 
the distance. Where independent distance measurements through 
measurements of the luminosity of the tip of the red giant branch 
are available, we adopt the distances and uncertainties listed in 
the source paper. 



4.1. Sub-samples 

As will be discussed extensively in Sect. I5.41 we will split our 
sample into two sub-samples, namely the "profile width" sub- 
sample and the "rotation curve" sub-sample. The profile width 
sub-sample contains the galaxies for which it proved impossi- 
ble to derive a maximum rotation velocity from either the tilted- 
ring analysis (V v f) or from the position-velocity diagram (V pv ). 
For this sample, therefore, only W20 and W50 measurements are 
available. Their properties are summarized in Table [2] The ro- 
tation curve sub-sample is summarized in Table [3] It contains 
the galaxies for which we were able to derive V v f and/or V pv (in 
addition to W20 and W50). 

In the Appendix, we show for all galaxies presented in this 
paper a summary panel consisting of moment maps, major and 
minor axis position-velocity diagrams and the global H 1 profile. 
We additionally show channel maps for the galaxies of the rota- 
tion curve sub-sample. 
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Table 2. Derived parameters for the profile width sub-sample 
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0.44 
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28 


19 


16 


10 


50 


D640-13 


10 56 13.9 


+ 12 0041 


958 


13.4 


-14.36 


0.55 


4.53 


38 


25 


30 


18 


48" 


D646-7 


12 58 40.4 


+ 14 13 03 


233 


2.1* 


-12.52 


0.84 


0.37 


35 


23 


26 


16 


55° 



optical /-band inclination from lPildis et al.1 i 19971) 

* distance as given in Karachents ev et all (2003) 
Notes: (1): galaxy identifier; (2, 3): central position taken from NED; (4): systemic velocity derived from the central velocity of the W50 profile; 
(5): distance in Mpc. If no referenc e is given, the dista nce is based on V sys (column 4) and a Hubble flow using H = 75 km s -1 Mpc -1 , including 
a correction for Virgocentric infall (Mould et alJl200CO ; (6): absolute V-band magnitude as given in lPildis et al .1 i 19971) . corrected to our distance 
estimates; (7): V-I color from lPildis et all {19971); (8): total Hi mass (in units of 10 7 Mq); (9): uncorrected width of the Hi profile at the 20 percent 
level of the maximum intensity. The values listed here for W20 are derived using the Hanning smoothed data; (10): as column 9, but for the IV50 
profile; (11): width of the W20 profile, corrected for the finite velocity resolution, and turbulent motion of the Hi gas; (12): as column 11, but for 
the W50 profile; (13): adopted inclination angle. The velocity widths in Cols. (9) - (12) are not corrected for inclination effects. 

Table 3. Derived parameters for the rotation curve sub-sample 
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54 
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63 ± 7 


13 
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12 52 21.9 


+21 37 46 
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-15.15 
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31.51 
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57 ± 5 


220 
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D631-7* 


07 57 01.8 


+ 14 23 27 


311 


5.5 4 


-14.50 


0.55 


14.68 


113 


89 


100 


78 


59 ±3 
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58 


53 



optical /-band inclination from lPildis et al.1 j 19971) 

* distance as given in Karachentsev et al. ( 2003) 
Notes: Cols. (1) to (12) as in Table[2] except that for the starred galaxies, the center position and V sys were derived kinematically; (13) inclination 
angle (optical inclination for D500-3, kinematic inclination otherwise). The uncertainties given for the kinematic inclination angles represent the 
scatter of the values of the individual tilted-rings. For D500-3, the average uncertainty of the kinematic inclination estimates is given; (14) average 
position angle of the rotation curve analysis, measured counter-clockwise from north to the receding side of the galaxy; (15) maximum rotation 
velocity from the tilted-ring analysis of the velocity field, corrected for inclination; (16) maximum rotation velocity from the position-velocity 
diagram, corrected for inclination. 



4.2. D500-2 (Data presented in Figs.U] \AJ\ and\AJ3 

D500-2 is also known as UGC 5716 and is classified as an Sm 
galaxy. We assume a distance of 17.9 + 3.4 Mpc. The global Hi 
profile of D500-2 (cf. Fig. I A. Il l shows the double-horned pro- 
file typical for spiral galaxies, and its velocity field indicates that 
D500-2 is clearly dominated by rotation. We fit ellipses using the 
zeroth moment map at a few representative intensity levels. The 
resulting inclination of the H 1 disk, corrected for the beam size, 
is /hi ~ 57°. The parameters from these isophote fits are used 
as initial estimates for the tilted-ring analysis. The width of the 
rings is set to 11". The systemic velocity as derived in the first 
tilted-ring fit shows only small radial variation. Its mean value 
is V sys = 1259 kms -1 , identical to the the central velocity of the 
W50 profile (1259 kms -1 ). Fixing the systemic velocity, we de- 
rive the dynamical center in a new fit by averaging the (xo, yo) 
values over the entire radial range. The resulting central position 
(cf. Table [3j is in good agreement with the optical center from 
NED and with that derived from our ellipse fits to the zeroth mo- 
ment map. Keeping the center position fixed for the subsequent 
tilted-ring fits, we derive the PA and inclination angle. The po- 
sition angle shows a linear decrease from 350° in the center to 
340° in the outer parts. 

The inclination is then determined in an additional fit by av- 
eraging the inclination values for r > 35". The resulting kine- 



matic inclination is /km ~ 57° + 6°, which agrees with the (beam 
corrected) inclination of the H 1 disk. In a last tilted-ring fit, we 
derive the rotation velocity by keeping all parameters except V 10t 
fixed to their previously derived values. After a gentle inner rise, 
the rotation curve of D500-2 reaches a flat part at a velocity of 
V v f ~ 68 kms -1 . The rotation curve, as well as the radial dis- 
tributions for PA and /km are shown in Fig. \T\ The maximum 
rotation velocity from the V-diagram is V pv ~ 68 km s -1 . 

4.3. D500-3 (Data presented in Fias.\A2}and \A.13[ 

D500-3 is a dwarf irregular (dl) at a distance of 18.5 + 3.2 
Mpc. Its velocity field (cf. Fig. IA.21 > indicates solid-body rota- 
tion throughout the entire disk, as confirmed by the V-diagram. 
We fit isophotes to different intensity levels in the zeroth mo- 
ment map and find an inclination of the Hi disk i m ~ 42°. This 
value is somewhat sm aller than the opti cal /-band inclination of 
/ opt ~ 55° as given in iPildis et al.l dl997l) . Because of the large 
H 1 beam size, we consider the optical inclination superior to the 
Hi inclination and use it in our further analysis. The solid-body 
rotation and the limited spatial resolution make it impossible to 
derive a well-determined rotation curve using tilted-ring mod- 
els. We are, however, able to get an estimate for the maximum 
rotation velocity by using the position-velocity diagram and find 
Vn V ~ 45 kms -1 . 
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Fig. 1. Tilted-ring analysis of D500-2. From top to bottom, the 
radial distributions of the rotation velocity, the inclination angle, 
and the position angle are shown. The data points in the center 
and bottom panel indicate the values derived for i and PA as free 
parameters. The solid lines indicate the values adopted to derive 
the rotation curve shown in the top panel. 

4.4. D512-2 (Data presented in Figs.\2\\KE and\AJ3 

D512-2 is a galaxy of Hubble type Sm. We assume a distance 
of 14.1 + 2.2 Mpc. The velocity field of D512-2 shows signs of 
solid-body rotation in its inner parts. We fit ellipses at a number 
of representative intensity levels in the zeroth moment map and 
derive an inclination of the H i disk of z'm ~ 48° (corrected for 
the beam size). Because the apparent size of D512-2 is small, 
we do not fit the center position or the systemic velocity using 
tilted-ring fits. For the center position, we adopt the estimate de- 
rived from the ellipse fitting, which is in good agreement with 
the position of the optical center. For the systemic velocity, we 
use the central velocity of the W50 profile. Keeping the center 
and the systemic velocity fixed, we derive an average position 
angle PA ~ 40° from the tilted ring fit, choosing the width of the 
tilted-rings to be 13". The position angle is then kept fixed for 
a subsequent tilted-ring fit to estimate i. Averaging the resulting 
inclination values yields z^n ~ 56° + 10°, whichis consistent with 
the values discussed above. Using the kinematic inclination, we 
determine the rotation curve in a last tilted-ring run leaving only 
V rot unconstrained. The rotation curve (cf. Fig. [2]) confirms what 
was already suggested by the velocity field: a solid-body rotation 
in the inner parts and a flat part in the outer regions. The rotation 
velocity of the flat part of the rotation curve is V v f ~ 35kms _1 . 
From the /jV-diagram, we estimate the (inclination corrected) 
maximum rotation velocity to be V pv ~ 37 km s _I . 

4.5. D564-8 (Data presented in Figs. and\Al3 

For the dwarf irreg ular D564-8 (also known as F564-V3, see 
de Bl oketal.fl 996). we assume a distance of 6.5+2.6 Mpc. The 
global Hi profile of D564-8 is asymmetric with its peak flux 
density towards the receding side of the galaxy. The asymmetry 
of the global H 1 profile can be traced also in the /? V-diagram. 
The maximum rotation velocity, obtained from the /^V-diagram 
is Vpv ~ 29kms~'. We use ellipse fitting to the zeroth moment 
map to derive initial estimates for our tilted-ring fits. We set the 
width of the tilted-rings to 12". In the first tilted-ring fit, we 
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Fig. 2. Tilted-ring analysis of D5 12-2. The layout is identical to 
that of Fig. □ 



derive a systemic velocity of V sys « 478 km s , which agrees 
well with the central velocity of the W50 profile. In a subsequent 
fit, we derive the dynamical center by averaging the values for 
Xq and yo over the entire radial range. The resulting center (cf. 
Tabled agrees to within a few arcseconds with the optical center 
from NED. In two subsequent tilted-ring fits, we derive the PA 
and inclination by averaging the individual measurements over 
the radial range, excluding the innermost (deviant) data point. 
The resulting inclination (z'ki n ~ 63° + 7°) is larger than both the 
one derived from t he H 1 disk (z'hi ~ 50°), and the optical incli- 
nation (/ pt ~ 35°, Pildis et al. ll997l) . However, an inclination of 
50°, or even 35° can be ruled out by our kinematic data. Keeping 
all parameters except the rotation velocity fixed to their best es- 
timates, we derive the rotation curve of D564-8. The maximum 
rotation velocity is V v f ~ 25kms~ 1 (see Fig. [3}. 
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4.6. D565-5 

The WSRT observations of D565-5 do not detect any emission 
at the reported position and redshift. We do, however, see strong 
emission from the nearby galaxy NGC 2903 at the edge of the 
primary beam. The agre ement of the apparent v elocity of this 
emission with that of the Schombert et al. ( 1997) detection sug- 
gests that the latter detection was simply NGC 2903 emission 
that was picked up with the larger single-dish beam. 

4.7. D572-5 (Data presented in F/a.lAB 

D572-5 is an irregular galaxy at a distance of 14.6+3.2 Mpc. Our 
H i data of this galaxy are barely resolved (see, e.g., the moment 
maps in Fig. IA.5b . which makes it difficult to derive a reliable H i 
inclination. We therefore adopt the value of the /-band inclina- 
tion / opt ~ 50° as given by Pi ldis et ail (I 19971) in further analysis. 
Due to the low resolution, substantial off-axis H i emission enters 
the beam, causing the large velocity range in the minor axis pV- 
diagram. Given these problems, a rotation curve or a V pv from 
the major axis position-velocity diagram cannot be reliably de- 
rived. The only usable indicators for the maximum rotation ve- 
locity are W20 and W50. 



4.8. D575-1 (Data presented in Fig.\A6} 

D575-1, also known as IC 3810, is classified as an Sm/Irr galaxy. 
We assume a distance of 10.0 + 3.6 Mpc. Its global Hi profile 
is well described by a Gaussian. The first-moment map and the 
position-velocity diagrams show only a small velocity range and 
it proved impossible to estimate either V v f or V pv . The inclination 
of the H 1 disk (derived using ellipse fitting and corrected for the 
beam size) is /hi ~ 53° and in rea sonable agreemen t with the 
optical inclination (/ opt ~ 61°) from lPildis et alj d 1997b (though 
the observed small velocity range in the velocity field and pV 
diagram casts some doubt on the validity of these values). 



4.9. D575-2 (Data presented in Figs. @|A3 and\AjE 

D575-2, or UGC 801 1, is a galaxy of the Hubble type Im, with 
an assumed distance of 12.2 + 4.5 Mpc. It was not possible 
to estimate V pv from the position-velocity diagram. The galaxy 
is kinematically lopsided, which can be seen in the differences 
in the velocity contours between the approaching and receding 
side. Using ellipse fitting, we determine the inclination of the H 1 
disk to be /m ~ 52 °, which is somewhat less inclined than what 
iPildis et all (fl997l) found optically (Z opt ~ 63°). With the Hi el- 
lipse fitting results as initial estimates for a tilted-ring analysis 
and adopting a width of the tilted-rings of 14", we determine a 
systemic velocity of (V sys ~ 774 km s -1 ). This agrees to within 
2kms _1 with the center of the W50 profile. Fixing the systemic 
velocity, we determine a position of the dynamical center (see 
Table O in excellent agreement with the optical center as given 
in NED. Leaving PA and / unconstrained, we find a gradual de- 
crease of the position angle from ~ 230° in the inner parts to 
~ 210° in the outskirts of D575-2. Fixing the PA to these val- 
ues, we derive the inclination by averaging over all data points 
with r > 0". The resulting inclination (hm~ 57° ± 5°) is halfway 
between the optical inclination and the one from the Hi disk. 
Keeping the inclination fixed to the kinematic estimate, we de- 
rive the rotation curve of D575-2 (see Fig. @}. After a linear 
increase in the inner parts of the galaxy, the rotation velocity 
reaches a flat part at ~ 74 km s _1 . 
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Fig. 4. Tilted-ring analysis of D575-2. The layout is identical to 
that of Fig. □ 



4.10. D575-5 (Data presented in Fig. 

D575-5 is classified as a dwarf irregular at an assumed distance 
of 7.7 + 3.1 Mpc. Its global Hi profile has a Gaussian shape 
and its W50 is the smallest of our s ample (W50 ~ 19k ms 1 for 
the Hanning smoothed data cube). iPildis et ail d 19971) estimate 
the optical inclination to be / opt ~ 66°, which is higher than the 
(beam corrected) inclination of the H 1 disk, for which we derive 
im ~ 50°. However, neither the velocity field, nor the major axis 
pV-diagram show clear signs of rotation. We were unable to de- 
rive a maximum rotation velocity from either the velocity field, 
or from the major-axis position-velocity diagram. The second- 
moment map shows values of 6-8 kms -1 throughout the entire 
disk. This suggests D575-5 is close to face-on and that a rela- 
tively large fraction of the profile width is caused by turbulence. 
For the inclination correction of W20 and W50, we use the in- 
clination of the Hi disk, as the appearance of D575-5 rules out 
the optical inclination of / opt ~ 66°, though we note that the ob- 
served small velocity range indicates that even this H 1 value is 
likely to be an overestimate of the true inclination. 

4.11. D631-7 (Data presented in Figs. El and \AJ7[ 

D631-7 i s also known as UGC 411 5 and is classified as a dwarf 
irregular. Karachen tsev et all d2003l) estimate its distance to be 
5.5 ± 0.6 Mpc using the luminosity of the tip of the red giant 
branch. The global Hi profile of D631-7 is single-peaked. Its 
velocity field is well-resolved and shows clear signs of rotation. 
The major axis pV-diagram shows indications of a flat rotation 
curve in the outer parts of the galaxy. Ellipse fitting yields an 
inclination of the Hi disk /hi ~ 57°. The ellipse fit results are 
used as initial estimates for the tilted ring analysis. The width 
of the rings is set to 12". In the first tilted-ring fit, we derive a 
systemic velocity of V sys ~ 311 kms -1 , which agrees well with 
the center of the W50 profile. The dynamical center is fitted in the 
second run with rotcur. The resulting center position is given 
in Table [3] and agrees to within 1" with the optical center as 
given by NED. The PA of D631-7 is obtained in a subsequent 
fit with rotcur and shows a gradual increase from ~ 318° in the 
inner parts to ~ 330° in the outer parts. Fixing the PA to these 
values, we derive the kinematic inclination of D63 1 -7 to be ~ 
59° + 3° by averaging the individual tilted-ring values for 50" < 
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r < 150". This is in good agreement with the inclination of the 
H i disk (/hi ~ 57°), b ut more face-on th an the optical inclination 
('opt ~ 66°, derived bv lPildis et al.l 19 971). Keeping all parameters 
except vrot fixed, we determine the rotation curve of D631-7 
(see Fig. O. It shows the typical solid-body rotation in the inner 
parts, but reaches a flat part at V v f ~ 58kms _1 . Assuming the 
kinematical inclination we derive V pv ~ 53kms _1 from the pV 
diagram. 
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Fig. 5. Tilted-ring analysis of D631-7. The layout is identical to 
that of Fig. Q] 



4.12. D640-13 (Data presented in Fig. \AM 

D640-13 is an Im/Sm type galaxy with an adopted distance of 
13.4 + 3.4 Mpc. Its global Hi profile has a Gaussian shape. 
Comparing the size and the shape of the beam with that of the 
moment maps shows that the galaxy is barely resolved. This 
heavily affects the derived inclination of the Hi disk (/hi ~ 71°) 
which is mu ch higher than the optical inclination (/ opt ~ 48°) by 
iPildis et all dl997j) . We therefore use the optical inclination for 
our further analysis. As neither the position-velocity diagrams, 
nor the first-moment map show clear signs of rotation, only the 
widths of the W20 and W50 profiles remain as a proxy for the 
maximum rotation velocity, which we correct using the optical 
inclination. 



4.13. D646-7 (Data presented in Fig. lATfl ) 

D646-7 ( orUGC 8091) is cla s sified as an irregular dwarf galaxy 
for which [Karachentsev et alj d2003l) estimate a distance of 2.1 + 
0.6 Mpc, using the tip of the red giant branch. The global H 1 
profile of D646-7 has a clear Gaussian shape. The galaxy lacks 
clear signs of rotation, both in the velocity field and in the major- 
axis position-velocity diagram. Thus, only W20 and W50 remain 
as indicators for the maximum rotation velocity. T he optical in- 
clination of D646-7, derived bv lPildis et all (Il997l) . is / opt -55°, 
which is in good agreement with the (uncorrected) inclination of 
the Hi disk (/hi ~ 51°). However, given that the galaxy is not 
well resolved, we use the optical inclination in our analysis. 
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Fig. 6. Stellar /-band mass-to-light ratio ( Y„) vs. (V - 1) color for 
two different popul ation synthesis models . The circles are based 
on the models from lBell & de Jong ( 2001 ) w ith a s caled Salpeter 
IMF, the triangles are based on the lPortinari et alj d2004l) models 
with a Kroupa IMF. Note that D631-7 and D640-13 have the 
same color: for clarity reasons we only plot one point. 



5. Analysis 

In the previous section, we have derived several estimates for 
the maximum rotation velocities of the galaxies in our sample. 
Before using them to construct a BTF, we discuss our choices 
for the stellar mass-to-light ratios (Y»), describe the corrections 
applied to the line width measurements and estimate the uncer- 
tainties in the rotation velocities and baryonic masses. 



5. 1. The choice of the stellar mass-to-light ratio T. 

One of the largest contributors to the scatter in baryonic mass 
is the assumed value for the stellar mass-to-light ratio Y„ and 
its uncertainty. Fortunately, the effect of Y* becomes less impor- 
tant when dealing with low-mass galaxies, as their stellar mass 
generally contributes less to the total baryonic mass than the gas 
mass does (for reasonable values of Y«). The Y* values adopted 
here are derived on the basis of two different population syn- 
thesis models using two diffe rent Initial Mas s Functions (IMFs). 
The first model is from Bell & de Jong (2001) and uses a scaled 
Salpeter IMF (Salpeter 1955). The second p opulation synt hesis 
model discussed here uses a Kroupa IMF dKroupal ll998) and 
is taken from lPortinari et all d2004l) . Figure [6] shows the derived 
stellar (/-band) mass-to-light ratios of our galaxies for the two 
population models as a function of color. 

It is immediately apparent that the Y„ from lBell & de Jong 

(2001) are consistently larger than the ones from lPortinari et al] 
d2004l) . Although this difference can amount to a factor close to 
three, the impact on the baryonic mass (i.e., the sum of the stellar 
mass and the gas mass) is rather small, as can be seen in Fig. [7] 
where, for each galaxy, we compare the ratio of the two different 
stellar mass assumptions against the baryonic mass implied by 
the respective Y» values. In the following analysis we will adopt 
the average Y» values implied by the two models. The individual 
mass-to-light ratios are listed in Table|4] 
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Table 4. Stellar, gas, and baryonic masses. 



i\J 


1 ^JJCll ) 


i yr oriinaii ) 




A/I 


AA f"R#ain 




/ AA \ 
V^baiV 










1 n8 A/i 
1U Mq 


1U Mq 


1 n8 A/* 
1U M G 


1U Mq 


(i) 


(2) 


(i) 


(4) 


(5) 


(6) 


(i) 


(8) 


D500-2 
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D500-3 


0.26 
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D575-1 
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0.28 


0.49 ± 0.21 


0.53 
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0.10 
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D575-2 


1.07 


0.45 


0.76 ± 0.31 


4.41 


1.02 


0.43 


5.13 


D575-5 


0.36 
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0.25 ±0.12 


0.60 
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D631-7 
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0.48 


0.18 


0.33 ±0.15 
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D646-7 


0.99 


0.41 


0.70 ± 0.29 


0.05 


0.09 
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N otes: (1): galaxy ide ntifier; (2): stellar /-band mass-to-light ratio based on Bell & de Jong (2001); (3): stellar /-band mass-to-light ratio based 
onlPortinari et al. (2004); (4): average of /-band T, values listed in columns (2) and (3). For the uncertainty in T„ we assume half the difference 
between the mass-to-light ratios of the two population models; (5): the gas mass in 10 s M Q ; (6): the stellar mass based on T» in column (2) in units 
of 10 s M Q ; (7): like (6), but based on T, in column (3); (8): the baryonic mass used in further analysis, derived from the gas mass in column (5) 
and the average of the stellar masses in columns (6) and (7). 
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(200lj)and|Portinari et al.l d2004l) vs. the implied baryoni c mass. 
Baryonic masses derived using the lBell&deJongl lfeOOl) stellar 
masses are shown as f i lled ci rcled; baryonic masses derived us- 
ing the iPortinari et al.l (12004) models as open circles. It is clear 
the choice of Y* has only a small effect on the baryonic masses. 

5.2. Line width corrections 

The rotation velocities derived using the profile widths W20 and 
W50, still need to be corrected for instrumental velocity resolu- 
tion and turbulence. T o correct for the instrumental resolution, 
we use the approach of lVerheijenl (Il997l) : 



W x = W° hs - C 



11 + 



(—)- 
\23.5) 



(5) 



where the subscript x refers to the chosen profile width measure 
(i.e., x — 20 for W20 and x — 50 for W50), R is the instrumental 
resolution in km s~' (cf. col. (9) in TableHJ, and C x is a constant, 
equal to C20 = 35.8 for the W20 profile and C50 = 23.5 for the 
W50 profile. The 'obs' superscript in W° bs denotes the observed 
profile width (given in cols. (9) and (10) of Tables|2]and[3]). 



In addition to broadening due to finite velocity resolution, 
we correct the velocity widths for broadening due to turb ulent 
motion of the neutral gas. Following T ully & Fouqud(ll985h . we 
use 



W 2 u = w 2 + w 2 

yr jc.turb " x T " t,x 











l-2e (^) 


-2W x W t>x 


1 -e 


( Wc,, ) 



(6) 



The subscript x again refers to the profile width measured at 
20 or 50 percent of the peak flux; W x indicates velocity widths 
already corrected for instrumental resolution following Eq. [5] 
The W CjA values represent the typical velocity widths where 
the shapes of the velocity profiles change from boxy (double- 
horned) to Gaussian; the constants W t , x indicate the amount 
by which turbulent motion of the neutral gas broad ens the Hi 
irofile . For the choice of W c , x and Wt, x , we follow IVerheijenl 
19971) . who assumes a turbulent motion of the gas with a ve- 
locity dispersion of 10 kms~' and values of W,2o = 22 kms~', 
W,,50 = 5 kms^ 1 W c .2o = 120 km s" 1 and W c , 5 o = 100 kms" 1 
(see Chapter 5 of I Verheii enll 1 9971 for an extensive discussion of 
these correction terms). The resulting corrected profile widths 
are given in Cols. (11) and (12) of Tables [2] and [3] Note that 
the profile widths listed in Tables [2] and [3] are not corrected for 
inclination. 



5.3. Sources of uncertainties 
5.3.1. Uncertainty in V max 

The uncertainty in V mdx depends on the measurement error in 
the observed velocity and that in the inclination value. For the 
velocity uncertainty we assume that we can measure the velocity 
with an accuracy of one channel width (~ 4kms _1 ). In addition, 
we assume that noise in the data or intrinsic variance within the 
galaxies introduce another channel width worth of uncertainty. 
The total uncertainty we assume is the quadratic sum of these 
two terms or 6 km s~' . 

For the uncertainty in the derived inclination angle, we dis- 
tinguish between galaxies with a kinematically derived inclina- 
tion, and those where we had to use the optical inclination or 
that of the H 1 disk. For the former category, we adopt the scatter 
of the inclination values of the individual tilted-ring models (see 
middle panels in Figs. [T]{5j as the inclination uncertainty (see 
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Col. (13) of Table [3). The average value of the uncertainties in 
these kinematically derived inclinations is 6° and is used as the 
uncertainty in the inclination angles of the remaining galaxies. 
The total uncertainty of the (inclination-corrected) maximum ro- 
tation velocity is then calculated using Gaussian error propaga- 
tion. 

5.3.2. Uncertainty of the baryonic mass 

The baryonic mass is the sum of the stellar mass and the gas 
mass: 

M h - dr = M stars +M gas . (7) 
The stellar mass is calculated through 

M stars = TilO- a4[ ^-^- 4 - 02] , (8) 

where T{ is the stellar mass-to-light ratio in the /-band, Mj = 
My — (V — I) is the absolute /-band magnitude, and 4.02 is the 
magnitude of the Sun in the /-band. The absolute magnitude of 
a galaxy depends through the distance modulus on its appar- 
ent magnitude and its distance. Thus, the stellar mass depends 
on three quantities: stellar mass-to-light ratio, distance, and ap- 
parent magnitude. For the uncertainty in the mass-to-light ratio, 
we use the values derived in Section [5TTI For the uncertainty in 
the distances, we use the values listed in Sect. [4] The apparent 
magnitudes of the galaxies in our sample were determined by 
Pildi s et all (1 19971) . who report a photometric accuracy of 0.05 
mag. This is insignificant compared to the influence of the un- 
certainties in the distance and in the stellar T„ and we therefore 
ignore the uncertainties of my for the uncertainty of the stellar 
mass. 

The second term contributing to the baryonic disk mass is 
M gas , the mass of the gas, which is given by: 

M gas = IAM H1 = 1.4 ■ 2.36 • 10 s D 2 j Sdv, (9) 

where Mhi is the total H i mass, D is the distance in Mpc, S is 
the total flux in mJybeairT 1 , and dv is the velocity resolution 
in kms -1 . The constant factor 1.4 corrects the Hi mass for the 
presence of helium and metals. Note that since H2/H1, the ratio 
between the molecular and the neutral hydrogen, is much lower 
in dwarf galaxies comp ared to luminous spirals dTavlor et al.l 
1998; iLerov et all 12005 ). we do not apply correction terms to 
account for molecular hydrogen. The uncertainty of the H 1 mass 
depends on the uncertainty of the total flux, and quadratically on 
the distance uncertainty. We focus again on the uncertainty in 
the distance, which is the dominant source of uncertainty here. 
Inserting Eqs. [8] and [9] into Eq.[7] we determine the uncertainty 
of the baryonic mass by assuming a Gaussian error propagation 
of the individual uncertainties discussed above. The stellar, gas, 
and baryonic masses are also listed in Table [4] 

5.4. The baryonic Tully-Fisher relation 

In this section we present the baryonic Tully-Fisher relation for 
the galaxies of our sample, using the different estimat es for V max 
derive d previously. As a reference, we use the work by McGaugh 
(2005). It presents the BTF for galaxies with well-determined 
rotation velocities between 50 kms -1 and 300 kms _I , and tests 
several methods to determine stellar mass-to-light ratios. The 
one yielding the smallest scatter in the BTF is ba sed on the mass- 
discrepancy-acceleration relation (MDAcc, see McGaugh 2004) 
and gives a BTF relation 

M bar = 50< ax , (10) 



with the baryonic mass in Mq and the maximum rotation veloc- 
ity in kms~'. 

For the construction of a BTF relation using our galaxies, 
we use the four different estimates for V max derived earlier. To 
recapitulate, these are: 

(a) V v f, obtained from a tilted-ring analysis of the velocity 
fields, 

(b) V pv , obtained from fitting the outer parts in major-axis 
position-velocity diagrams, 

(c) Vwio, obtained from ^W2o, the rotation velocity as de- 
rived from half the (corrected) width of the H 1 profile at the 
20 percent level of the maximum flux, 

(d) Vwso, obtained from hWso, the rotation velocity as de- 
rived from half the (corrected) width of the H 1 profile at the 
50 percent level of the maximum flux, 

These velocities are all corrected for inclination. Additional cor- 
rections have been applied to VV20 and Vwso (as described in 
Sectionl5T2l). 

In Fig [8] we show the baryonic Tully-Fisher relation for the 
galaxies of our sample using the different velocity estimates. 
The uncertainties are discussed in Section 1531 In Fig [8] we dis- 
tinguish between galaxies in which clear signatures of rotation 
were detected in the position-velocity diagram and/or the ve- 
locity field (y pv and/or V v f available in addition to VV20 and 
Vwso', the rotation curve sub-sample), and galaxies for which 
no clear rotation was detectable (only Vwio and Vwso available; 
the profile-width sub-sample). For all four measures of V max , the 
galaxies in the rota tion curve sub-sampl e are consistent wi th the 
BTF as derived in lMcGaughld2005l) and lStark et ail (120091) . 

A few galaxies in the profile-width sub-sample are also con- 
sistent with a line-width-based BTF, but the majority of these 
galaxies are found to the left of the extrapolated BTF. This may 
indicate an increased scatter at the low line-width end of a line- 
width based BTF. Alternatively, the galaxies may be more face- 
on than indicated by the optical or H 1 inclinations (e.g., D575-5), 
or they may not be supported by rotation. Without independent 
and/or resolved measures of the rotation velocity and the kine- 
matic inclinations, it is difficult to say anything further on these 
galaxies. 

In the following we therefore restrict our analysis to those 
galaxies which have a well resolved velocity field which allowed 
us to derive a rotation curve or — at the very least — a maximum 
rotation velocity using the posit ion-velocity di agr am (consis- 
tent w ith the procedure used in lMcGaughll2005l and IStark et ail 
2009). We emphasize that this is our only selection criterion. We 
do not preferentially select against galaxies that are not on the 
BTF, (in fact, one of the profile-width-only galaxies (D572-5) is 
right on the BTF), and we include all galaxies for which V pv or 
V v f could be determined regardless of their position in the BTF 
diagram. 

5.5. The scatter of the BTF 

The panels in Fig. [8]clearly show that the data points from our 
rotation curve sample agree well with the BTF relation from 
iMcGaughl (120051) as given in Eq. \W\ This relation was derived 
using resolved observations of the flat parts of the rotation curves 
of the sample galaxies. Including our five galaxies with mea- 
sured values of V w { yields a revised BTF of the form 

M bar = 58V^ X 7 . (11) 

This is remarkab ly close to the rev ised BTF presen ted in 
IStark et al.l d2009l) . which extends the iMcGaug h (2005) BTF 





V W20 ( km S ) 




10 1 1CT 
V w , n (km s" 1 ) 



Fig. 8. The baryonic Tully-Fisher relation using four different estimates for V m3X , namely V v f from the tilted-ring analysis of the 
velocity fields (top left), V pv from the major-axis position-velocity diagram (top right), Vwio, i-e-, half the W20 (bo ttom left), and 
Vw5Q, i.e., half the W50 (bottom right). As a refer ence, we show th e BTF relation of the analysis of dMcG augh 2005, dashed line). 
In the top panels, the star-dominated galaxies of McGaugh] d2005l) are shown as small square s and the gas-d ominated galaxies of 
Star k et aTT ( 2009b are shown as small circles. In the lower panels, the velocity widths from Bothu^ d 19851) (small squares) and 
McGaugh et aDd2000) (small circles) are shown for comparison. Our new data are show as large circles. Galaxies for which no 



rotation is discernable are shown as open circles (with line-widths only) and those with rotation det ected in the po sition-velocity 
diagram are shown as filled circles. The latter fall along the extrapolation of the BTF relation found by McGaughld2005l) while some 
of the former fall off the relation. These galaxies may not be supported by rotation, may have incorrect inclinations, or could simply 
reflect the increased scatter in the relation when line- width is used as a proxy for rotation velocity. 



with an independent sample of low-mass gas-rich galaxies. They 
find 

M bar = 62V^. (12) 

The three fits given in Eqs.[lO],[TT|and[T2]are statistically iden- 
tical. 

For comparison, Bell & de Jond d2001l) find a slope of 3.5 
for a subsample of the McGaugh (2005) sample, but choose dif- 
ferent Y» values than lMcGau gh (2005) does. Using identical Y» 



values yields i dentical slopes. Work b y lGeha et al.1 (|2()06) gives 
a slope of 3.7, |Pe Riicke et alj d2007l) find a slope of 3.15, but 
this in cludes early-type galaxies. An analysis by iMever et~ai1 
( 2008) of H i-selected galaxies gives a steeper slope of 4. 35 . The 
majority of these studies deal with massive galaxies, dominated 
by stars, and are therefore sensitive to assumptions on Y*. The 
new data presented here provide a test of the slope extrapolated 
from the fit to higher mass galaxies that do not overlap with this 



C. Trachternach et al.: Baryonic Tully-Fisher relation 



11 



samp le, and p r efer a slope of ~ 4. This is also the conclusion 
of IStark et al.l (120091) who calibrate the BTF using only gas- 
dominated galaxies, independent of assumptions on T». They 
find a slope of 3.94 (Eq. [T2"l>. If a single, unbroken BTF exists, 
this can only imply that the slope at the high mass end must also 
be close to 4. 

We now return to the conclusions of iFranx & de Zeeuwl 
(119921) . who use the scatter of the TF relation to put constraints 
on the ellipticity of dark matter halos. They argue that if the po- 
tential in the plane of the disk is elongated, then the different 
viewing angles will cause scatter in the TF relation. In return, 
the scatter of the (B)TF relation can be used to put an upper 
limit on the e llipticity of the potential i n the plane of the disk. 
According to IFranx & de Zeeuwl (|1992|) . a TF relation with a 
scatter of 0.31 mag (0.46 mag when photometric inclinations 
are used) constrains the ellipticities of the potentials to be be- 
low 0.10. As it is unlikely that all the scatter in the TF relation 
is due to different viewing angles, they argue that an elongation 
between and 0.06 is r easonable. 

Since the work of IFranx & de Zeeuwl (|1992|) . the quantity 
and quality of the data have improved significantly and we can 
now trace the BTF relation over a large range of galaxy masses 
and rotation velo cities. The scatter in the optimum BTF from 
iMcGaughl (120051) is cr = 0.098 dex, or 0.25 mag. Including 
our galaxies results in a slightly larger scatter (<x = 0.13 or 
cr M = 0.33 mag) — mainly because the distances to our galax- 
ies are more uncerta in. This scatter is only ma rginally larger than 
the scatter on which lFranx & de Zeeuwl d 19921) based their upper 
limits for the ellipticity of the potentials in the disk plane. 

Their arguments equally apply to our sample. As the small 
scatter in the BTF holds down to rotation velocities of ~ 
25 km s _1 , it follows that the elongation of the potentials of these 
ga laxies cannot be very large , and similar limits as those derived 
bv lFranx & de Zeeuwl (fl992). also apply to our galaxies. The up- 
per limit on elongation of the potential of 0.06 is consistent with 
the results derived by [Trachternach et all d2008l) who use har- 
monic decompositions of the velocity fields of 1 8 nearby well- 
resolved galaxies fro m the THINGS survey dWalter et al.ll2008t 
Ide Blok et al J 120081) to put constraints on the elongation of the 
potential in the plane of the disk. They find that the the aver- 
age elongation of the potential is small (0.017 + 0.020), particu- 
larly when compared to what is found in CDM simulations (cf . 
Fren k et dll98"alHavashi & Navarroll2006l) . lTrachternach et alj 
d2008l) also find that the elongation does not increase towards the 
center of the galaxies. The tightly correlated BTF relation pre- 
sented in this paper constrains the ellipticity of galaxy potentials 
using a different method, but reaches similar results. 

6. Conclusions 

We present new H i observations of a sample of low-mass dwarf 
galaxies and use these to explore and extend the baryonic Tully- 
Fisher (BTF) relation at low rotation velocities and galaxy 
masses. We present and discuss several estimates for V max , the 
maximum rotation velocity. For galaxies where a clear rota- 
tion signal could be detected, the different estimates are in good 
agreement and the BTF relations based on them are equally well- 
constrained. We discuss the choice of stellar mass-to-light ratio 
(T») and show that its choice is not crucial for the extreme dwarf 
galaxies in our sample, since for them, the stellar mass generally 
contributes less to the total baryonic mass than is the case for 
luminous high-mass galaxies. 

The small scatter in the BTF presented here (cr M = 
0.33 mag) puts strong constraints on the ellipticity of the po- 



tential in the plane of the disk o f the galaxies. Our result s 
are in agre e ment with t hose from IFranx & de Zeeuwl d 1992b . 
IFranx et aU d 19941) . and iTrachternach et all d2008l) . indicating 
that, at least in the plane of the disk, galaxy halos are not very 
elongated. The small scatter in the BTF over almost 5 orders of 
magnitude of baryonic mass indicates it is a fundamental rela- 
tion which tightly couples the visible baryonic matter and the 
dark matter content of galaxies. 
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Appendix A: Atlas 

The Appendix contains summary panels for all galaxies of our 
sample (Figs. lATTlA.l It and channel maps for the galaxies from 
the rotation curve sub-sample (Figs. |A~12llA. 171 ). 

The summary panels consist of two rows with three panels 
each, and contain the following maps: 

Top row: Left panel: zeroth moment map with grayscales. 
Unless mentioned otherwise, grayscales run from a column den- 
sity of n m = 1 • 10 19 cm" 2 (white) to n m = 2- 10 21 cm" 2 (black). 
The 3cr level is indicated by the black contour. Middle panel: 
first moment map. The systemic velocity is indicated by the thick 
contour, the contours are spaced by 10 kms" 1 . The approach- 
ing side can be identified by the light grayscales and black con- 
tours, and the receding side by dark grayscales and white con- 
tours. Right panel: second moment map. Grayscales run from 
2 to 40 kms" 1 . Unless mentioned otherwise, the contours lev- 
els are given at 5, 10, and 15 kms" 1 . For all moment maps, the 
beam size is indicated in the bottom right corner. Additionally, if 
a kinematic center has been independently derived, it is indicated 
in all moment maps by a cross. 

Bottom row: Left panel: Major axis position-velocity dia- 
gram. The position angle of the slice indicated in the top-left 
corner of the panel. Grayscales run from 2<x to 30<x, and con- 
tour levels are given at 2cr + n x Act (i.e., 2,6,10,14,18 ... <x). 
The dashed line indicates the systemic velocity. If a V pv has 
been derived, the resulting velocities at either side of the rota- 
tion curve are indicated by arrows. Middle panel: Minor axis 
position-velocity diagram. Grayscales and contours are identical 
to the major axis position-velocity diagram. Right panel: Global 
H i profile of the Hanning smoothed data cube. 
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Fig.A.l. Summary panel for D500-2. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3<x contour in the zeroth moment map corresponds to a flux of 5.09 mJy, which translates to a column density limit of n HI = 
8.8 • 10 19 cm 2 . 
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Fig.A.2. Summary panel for D500-3. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 4.51 mJy, which translates to a column density limit of nm = 
5.2 • 10 19 cm" 2 . 



14 



C. Trachternach et al.: Baryonic Tully-Fisher relation 



25" 



14 h 33 m 20 s 15 s 



25" 



R.A. (2000.0) 
14 h 33 m 20 s 



15 s 



25" 



14"33 m 20 s 15" 



50 

27° 0' 

26°59'10" 

26°58'20' 
900 

T 

V) 

E 850 
> 

800 




- PA = AOS, 





--"PA = 310° 






50 



27" 



- 26°59 10 



26°58 20 



0.03 



- 0.02 



- 0.01 



+ 1 -1 

Offset (arcmin) 



+ 1 -1 

Offset (arcmin) 



800 850 

V hel (km s- 1 ) 



900 



Fig.A.3. Summary panel for D512-2. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3<r contour in the zeroth moment map corresponds to a flux of 5.50 mJy, which translates to a column density limit of n HI = 
3.3 • 10 19 cm" 2 . 
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Fig.A.4. Summary panel for D564-8. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 4.99 mJy, which translates to a column density limit of nm = 
4.9 • 10 19 cm" 2 . 
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Fig.A.5. Summary panel for D572-5. A general description of the panels and levels is given at the beginning of the Appendix. 
Grayscales in the zeroth moment map run from a column density of n m = 1 • 10 cirT 2 (white) to n HI = 8 • 10 20 cm (black). 
The 3cr contour in the zeroth moment map corresponds to a flux of 12.47 mJy, which translates to a column density limit of 
n HI = 1.4- 10 19 cm- 2 . 
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Fig.A.6. Summary panel for D575-1. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3<x contour in the zeroth moment map corresponds to a flux of 9.81 mJy, which translates to a column density limit of n m = 
1.3 • 10 19 crrT 2 . The contours in the second moment map are at 7.5 kms -1 and 9 kms" 1 . 
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Fig.A.7. Summary panel for D575-2. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 5.58 mJy, which translates to a column density limit of n HI = 
2.0 • 10 19 cm 2 . 
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Fig.A.8. Summary panel for D575-5. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 8.82 mJy, which translates to a column density limit of nm = 
0.6 • 10 19 cirT 2 . The contours in the second moment map are given at 6 km s _1 and 7.5 km s _1 . 
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Fig.A.9. Summary panel for D631-7. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 11.06 mJy, which translates to a column density limit of 
n HI = 3.7 • 10 19 cm" 2 . 
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Fig. A.10. Summary panel for D640-13. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3<x contour in the zeroth moment map corresponds to a flux of 1 1 .49 mJy, which translates to a column density limit of 
n m = 3.3 • 10 19 crcT 2 . The contours in the second moment map are given at 7.5 kms~' and 10 kms -1 . 
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Fig. A.ll. Summary panel for D646-7. A general description of the panels and levels is given at the beginning of the Appendix. 
The 3cr contour in the zeroth moment map corresponds to a flux of 7.92 mJy, which translates to a column density limit of n m = 
0.7 • 10 19 cirT 2 . The contours in the second moment map are given at 7.5 km s _1 and 9 km s _1 . 
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Fig. A.12. Channel maps of D500-2. Grayscales run from -0.02 to 7 mJy. Every second channel is shown. 
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Fig. A.13. Channel maps of D500-3. Grayscales run from -0.02 to 5 mJy. Every channel is shown. 
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Fig. A.14. Channel maps of D512-2. Grayscales run from -0.02 to 8 mJy. Every second channel is shown. 
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Fig. A.15. Channel maps of D564-8. Grayscales run from -0.02 to 7 mJy. Every channel is shown. 
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Fig. A.16. Channel maps of D575-2. Grayscales run from -0.02 to 10 mJy. Every fourth channel is shown. 
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Fig. A.17. Channel maps of D631-7. Grayscales run from -0.02 to 25 mJy. Every second channel is shown. 



